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One and Two Dimensional Spin Systems in the Regime Close to Deconfinement of 
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Based on the Majumdar-Ghosh chain we construct several spin models which allow us to inves- 
tigate spinon dynamics in the regime close to deconfinement of spinons. We consider the Ji — J2 — S 
model, two coupled Ji — J2 chains (ladder), and a 2D array of coupled Ji — J2 chains. Using the 
picture of two spinons interacting with a string confining potential we calculate the singlet-triplet 
splitting, magnetic structure factor, tunneling amplitude of two spinons and the excitation spectra 
for the ladder and the array. 
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(D 

^^ . Interest in low dimensional systems with dinierization has greatly inxireased injxecent years due to severaLexperimen- 

00 ■ tal findings. Compounds such as (VO)2P207y, Cu(N03)2 ■ 2.5H2OH, CuWOj and Cu2(C2Hi2N2)2Cl4y effectively 
behave as one-dimensional chains with intrinsic dinierization. Further motivation for the interest comes with the 
discovery of a spin-Peierls transition in the dimerized, weakly coupled, ladder system, CuGeOsH which now has a 
well determined magnetic structured. Finally,_the recent discovery that the ladder system, Sro.4Ca13.gjCu24O41.g4, 
becomes a superconductor at 12 K and 3 GPaEl has renewed the idea, originally suggested by Andersora, that spin 
^^ '■ chain dynamics is linked with the mechanism of high-Tc superconductivity. 
04 ■ The interest in systems with dinierization is not restricted to the ID case. There are a number of 2D theoretical 

Cn ' models with induced, or spontaneous dimcrimtion, such as the two-layer Heisenberg model (see, e.g. Ref.EI), dimerized 
^^ , Heisenberg modelalj and the Ji — J2 modelllil. Some of the models can be relevant to real materialalj. 

A very important question relevant to any disordered quantum spin system is: what are the elementary excitations 
of the system? Until recently, common wisdom was that there are two possibilities: 1) elementary excitations have spin 
1 (elementary triplets) or 2) elementary excitations have spin 1/2 (spinons). See RefEj for a review. However, now 
it is becoming clear that the spectrum of excitations of disordered quantum spin systems is much more complicated. 
The spin 1/2 ID Heisenberg chain is a very important example of a pure spinon excitation spectrumtj. Some 
'^ ' generalizations of this model, such as the frustrated Heisenberg chain, also have pure spinon excitations. However, 
^ introduction of external dinierization changes-the spectrum drastically: the spinons become confined and the spectrum 

O consists of multiple singlet and triplet stateai3'E3. On the other hand, at very strong dimerizations, one would expect 
pure triplet excitations. lChia.is, however, not the case. It has been demonstrated that there are multiple bound states 
of the triplet excitationstlTtJ. Morepwer, the bound state can sometimes have energies very close to the "elementary 
triplet" , as in the 2D Ji — J2 modelESI. The complexity of the spectra of quantum spin systems is somewhat similar 
to the complexity of the meson spectra in Quantum Chromodynamics. In reality, the physics is also very similar. 
C^ ' The purpose of the present study is to investigate the excitation spectra of the spin systems which are close to 

deconfinement of spinons. As a building block we use the frustrated Ji — J2 spin 1/2 chain near the Majumdar-Ghosh 
point {J2/J1 — 0.5). Spinons are deconfined in this model. The models we consider have bound states of arbitrarily 
large size, meaning that they are close to deconfinement. The rest of the paper is organized as,follows. In Section II 
we consider a Ji ~ J2 — S chain at J2 ~ 0.5Ji. Following the approaches suggested by Afflecktj and Uhrig et alli3 we 
calculate singlet-triplet splittings and spin structure factors (quasiparticle residues for triplet excitations). In Section 
III we consider two coupled Ji — J2 chains (a ladder). The most interesting part here is the tunneling of the spinon 
from one leg of the ladder to another. In Section IV we consider two-dimensional array of the Ji — J2 chains. We 
calculate the confinement radius and demonstrate that the low energy dispersion of the triplet excitation is isotropic. 
Section V presents our conclusions. 

II. ONE DIMENSIONAL ,/i - .h - 5 MODEL 

The Hamiltonian for this model is of the form 

^ = -^ E [(1 + (-'^)''^') S.S.+i + aS,S,+2] (1) 
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where Si represents spin 1/2 at the site i. We set Ji = J, J2 = aJ and denote the degree of dinterization as S. 
Consider first the case of zero expHcit dimerization (5 = 0. It is well known that sd a < ac ~ 0.24lLil'E3 there is a 
unique ground state and that the excitation spectrum is gapless. On the other hand, at a > ckc, the ground state 
has spontaneous dimerization and is therefore doaJply degenerate, with a gap in the excitation spectrum. It has been 
shown by Majumdar, Ghosh and van den Brocket that, at a = 1/2, all quantum fluctuations are canceled out and 
the ground state is particularly simple: the product of exact spin dimers, 
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[j-2,j-l][j,j + l][j + 2,j + 3] 



(2) 



where we denote the singlet pair formed from adjacent sites, i and j, by [i,j]. The singlet chain shifted by unity with 
respect to Eq. (0) is also a valid ground state but the degeneracy is broken when 6^0. Due to the simplicity of the 
MG point we largely restrict our analysis in the present work to the case of a = ^ ■ 1—1 

Another important development was the analysis of the excitations on a MG chain by Shastry and SutherlandcJ. 
The simplest excitation consists of a pair of propagating topological kinks (spinous), see Fig. 1. 
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FIG. 1. Propagation of a pair of spinons (arrows mark the unpaired spins) on the MG chain. Ovals represent singlet pairs. 
Note that the spinons may not pass through each other due to the dimerized nature of the chain. 



Using the wave function corresponding to Fig. 1 as a variational function one can easily find the spinon dispersion 

(3) 
The energy gap is therefore A — J/8. The dispersion (0) can be expanded near minimum, q = k — tt/2, as 
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<^k = -^J + -Jcos2k. 
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with effective spinon niass-|rn = Q.b/ J. In spite of the variational nature of the dispersion, (||) agrees very well with 
exact numerical resultsOES. The scattering continuum for a pair of spinons with total momentum K is shown on Fig. 
2. 




FIG. 2. The scattering continuum for a pair of kinks with total momentum K as found in RefEj (hashed regions) and the 
dispersion for the triplet with unity spinon separation and a = 5 (dashed line). J has been set to unity. We see that the triplet 
configuration is a bound state at A' = ^. 



Two spinons with parallel spins interact with one another. There is a very simple way to find this interaction. Let 
us consider the state with spinons at nearest sites, see Fig. 3. 
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FIG. 3. Triplet hopping with unity kink separation. 

This corresponds to the triplet excitation of the spin-singlet dimer. The excitation energy of the triplet is J. In 
general the triplet can hop to the nearest bond, but the hopping matrix element 0.25J(1 — 2a) vanishes at the MG 
point. Thus, the energy of the triplet shown in Fig. 2 by the dashed line is independent of momentum: Etripiet = J ■ 
Certainly this is not an exact result. Quantum fluctuations ensure there are some corrections to this energy, but they 
are very small. The triplet does not exist as a stationary state when its energy is inside the two spinon continuum. 
However, we see from the Fig. 2 that there is a window where it does exist. This state has been found earlier in 
Ref.EJ using a more sophisticated method. The triplet exists because of the effective attraction between the spinons 
at nearest sites which we can write as 

Veff {Xi ,X2) = C5xi + l,x2 , (5) 

where xi and X2 are the spinous' coordinates which take integer values. To find the constant C let us write down 
wave function of the triplet with momentum K as 

^k{xi,X2) oc e*^^4,i ex e'^^ ^ e*M--i)^ (6) 

fe 

where X — (xi + X2)/'2 is the center of mass coordinate and a; = X2 — a;i is the relative coordinate of the two spinons. 
In the bound triplet, x = 1. Eq. (|6|) can be transformed to 

^Kixi,X2) (X ^ ^~^k^^(f+k)x,^^(f-k)x2_ (7) 

k 

In this wave function all momenta contribute with equal weight, therefore the average kinetic energy is 

-[e{K/2 + k) + e{K/2-k)] = p, (8) 

-7r/2 ■^ 4 

where e(fc) is the spinon dispersion (|3|). For two spinons at adjacent sites the average value of the effective interaction 
(|) is (*|Ve//|*) = C. However, we know that the total energy of triplet is J = {^\Ekin\'i') + (^IK//!*)- From here 
we find that C — — t^ and hence the effective spinon-spinon interaction in the triplet channel (S = 1) is 

Veff{xi,X2)^--JSx,+l^a:2- (9) 

No such potential exists for the singlet {S ~ 0) state of two spinons. 

Up till now we considered zero explicit dimerization, S — 0, which results in two degenerate ground states and 
deconfinement of spinons. Let us now consider very small but nonzero dimerization. The degeneracy between the 
ground states is removed and hence creation of two spinons creates a string with tension 3(5 J/4. This is confinement. 
For low energy excitations one can use the quadratic approximation (j^ and, hpn ce the wave function of the relative 
motion of two spinons {x — X2 — xi) obeys a Schrodinger equation, see Refs.tSEJ 

{E-2A)^{x)=- — -^^j{x) + ^x^P{x), (10) 

fj, is the reduced mass, /i — m/2. Boundary conditions are important here. The boundary condition at large x is 
obvious: ijj{oo) — 0. Boundary conditions for small x are different for triplets and singlets 

triplet : ^\-l) = 0, (11) 

singlet : ^'''(+1) = 0. 



The condition for the triplet reflects the fact that the spinons (kinks) cannot penetrate through each other, see Fig. 
1. This is also true for the singlet, but in this case there is an additional condition: excited states must be orthogonal 
to the ground state. Therefore, we assert that "0^(1) = since the singlet with unity kink separation is identical to 
the ground state. Note that the bound st|ai,eS| have very large sizes {x 3> 1) and , hence, in the first approximation, 
one can replace (O) by "0(0) = 0, see Rcfs.ll3'ta. We must use (O) since we intend to consider triplet-singlet splitting. 
Equation ( |lO| ) with boundary conditions (|ll| ) gives spectrum 



E(t-^^) = 2A - 2znJ i -S 
where z„ {n — 0, 1, 2...) are the zeros of the Airy functior 

zo = -2.388, zi = -4, 
The wave function has the form 
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(12) 
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where ^ = (3(5/8) ~^/^ is the typical size. The signs ± in Eqs. ([l^ ) and ( [l4| ) corresponds to triplet and singlet states. 
The upper sign always corresponds to the triplet. 

One may conclude from Eq. (12) that the singlet-triplet splitting is |(5J. However, this is not quite correct because 
Eq. (ITH) does not take account of the attraction (0) in the triplet channel. This attraction decreases the triplet energy 

by 



(V''|K//IV'*> ^ 

Together with (112) this gives the singlet-triplet splitting 
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which is independent of n. Certainly this is not an exact result because some small quantum fluctuations were 
n^ftljected in our derivation. Nevertheless it agrees very well with numerical data obtained recently by S0rensen et 
1 The data are presented in Table I. 
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Technique 



No. of sites 



0.050 
0.025 
0.005 



0.092 
0.049 
0.011 



(a) 
(a) 
(b) 



512 
512 

28 



TABLE I. The value of the singlet-triplet splitting obtained via (a) DMRG techniques described in Ref.ES and (b) exact 
diagonahzation (J has been set to unity). The finite size effects are negligible. Data provided by Erik S. S0rensen. 



Our result (^^ is valid for small momenta of the bound states, K Ki (or K r^ n). Nevertheless it practically 
coincides with|-splitting which is-known a,t K = tt/2. The singlet and triplet excited states at (5 = and K = 7r/2 are 
known exactlyc3. One can showES that these excited states remain eigenstates when we introduce explicit dimerization 
and that the triplet and singlet have energies (1 + 6)J and (1 + 36) J, and so the splitting is 26 J. 

In conclusion of this section let us the discuss magnetic structure factor or spectral weight of the triplet excitations. 
The external perturbation is 
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(17) 



where Sj^x is the spin raising-|Operator at the site x. Then the amplitude of the spin triplet creation at a spontaneously 
dimerized bond is (sin ■^) Eil. However this triplet is only a virtual state. To find spectral weight of the stationary 
state f/'n we have to project the virtual state onto the stationary one. This is proportional to the amplitude for two 
spinous to be at neighbouring sites. Hence 



Z„ = sm — h/'„(l) =7TSm — = -(5sm — . 



Our consideration of the bound states is valid only at small K or at K 
The spectrum of triplet excitations at i^T = tt is shown in Fig. 4. 



(18) 



TT and hence (|l8|) is valid only at these K. 
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FIG. 4. Schematic plot of the spectrum of excitations described by Eqs. (|12|),(|13|),(|18|) and (|l9|). We note the gap < u < 2A, 
the set of discrete bound states for 2A < u) < 4A with spectral weight indepenuent of n and tne continuum for ui > 4A. 
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Schematic plot of the spectrum of excitations described by Eqs. (|l2|) , (113) 



There are no excitations at w < 2A = 0.25J. At a; > 4A = 0.5J the triplet excitation can decay into pairs 
(triplet+triplet or singlet+triplet) and therefore the spectrum is continuous. At 2A < a; < 4A there are a number of 



discrete bound states. The spectral weight of any of these states is given by Eq. (|18|). The number of states, Umax, 
depends on 6, such that for a small 6 there is a large Umax- One can easily find Umax from Eq. (12). For example, 
a.t 6 > 0.14 there are only two states: n = and n = 1. At 6 = 0.01 the number is rimax = 6, and a.t 6 = 0.001 the 
number is Umax = 30. At extremely small 5 the asymptotic formula is valid: rimax ~ 0.025/^. Finally the splitting 
between the nearest peaks at n ::^ 1 according to (|2|) is equal to 
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1.95J(5^/3 

(n + 3/4)1/3 ■ 



(19) 
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The spectral weight in Eq. ( p^ is not quite consistent with that obtained using a series expansion methodcS 
possible reason for this is that convergence of the series may be very slow due to the dense spectrum (see Fig. (|l|)). 



III. TWO COUPLED MAJUMDAR-GHOSH CHAINS (THE LADDER) 



Let us consider two MG chains coupled to the ladder 

H = ,/^ {[S.S.+i + 0.5S,S,:+2] + [S:s^+i + 0.5S,:S,:+2] } + ^± E ^^^^ 



(20) 



with J± <C J. At J^ = there are four degenerate ground states. Two of them are topologically different and are 
shown in Figs. 5a and 5b. 
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FIG. 5. (a) Two chains in the 'normal' state and (b) the 'ahernate' state. 

These will be called "normal" and "alternate" states respectively. At Jj^ ^^ the degeneracy is lifted and the 
"normal" state corresponding to Fig. 5a becomes the true ground state. To see this, one needs to consider the second 
order (in J±) correction to the ground state energy. To calculate this correction it is convenient to use a localized 
triplet representation for the virtual excitation. Virtual admixtures of two triplet excitations to the "normal" state 
are shown in Fig. 6. 
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FIG. 6. Virtual excitations of the 'normal' ground state to the state with two triplet pairs. 

Keeping in mind that each localized triplet has energy J, we can easily find the correction AEnormai = — 3 J^/8JE3. 
Quite similarly the correction for the "alternate" state due to virtual excitations shown on Fig. 7 gives AEaiterate = 
-3Ji/I6J. 
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FIG. 7. Virtual excitations of the 'alternate' ground state to two different states, each with two triplet pairs. 
Thus, the "normal" state has a lower energy with the difference 

3J2 

loJ 

per two rungs along the ladder. We should note that this calculation is not exact. Because of quantum fluctuations 
in the virtual states there is a small correction to Eq. ( pf| ) which is also proportional to J^. Similar to the case of 
spinon dispersion (|^), we neglect this correction. 

Two spinon excitations on the one leg is shown in Fig. 8. 



FIG. 8. Two spinon excitations on one leg of the ladder. 



It is clear from ( pi| ) that the tension in the confining string is T = 3J^/32J and, therefore, the system is described 
by Schrodinger equation ( |lO|) with substitution 

3W 3 {JjAf 

4 2 J ' ^ ^ 

According to Eq. (q3) the wave function for two spinous with total momentum K (along the ladder) is of the form 

its), ^ e'^^ Ai{{x± l)/S. + z„) , , 

il^\i'{xi,X2)^—r=- r^ , (23) 



1/3 /^jx2/3 



where X — [xi + a;2)/2, x = a;2 — xi, and the index L = 1, 2 describes the leg on which the spinons are moving: L = 1 
corresponds to the upper leg and L — 2 corresponds to the lower one. Let us consider now the specific case of triplet 
excitation. There is a possibility for the triplet to tunnel from one leg to another. The mechanism of the tunneling is 
shown in Fig. 9: First the two spinons must come to adjacent sites, then the perturbation, J_lSS', swaps the triplet 
from one leg to another, and then the spinons propagate along the second leg. 
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FIG. 9. The mechanism for tunnelling of two spinons. The spinons on the top leg of the ladder must first come together 
before they can tunnel to the lower leg. After tunneling, the spinons are free to propagate along the lower leg. 

Simple calculation shows that the tunneling matrix element is equal to 

r = {i^l2\T.^.S'M,) = ^ =6^:^. (24) 

It is interesting that this is independent of the indexes n and m. According to (n3) (with substitution (p2)) the 
energy splittings between states with different n is ex J^ . The tunneling matrix element, r oc J^, is much smaller 
and therefore mixing within one n-level is independent of the others. Effectively, at each n we have a two degenerate 
level system with tunneling r between them. The stationary states have a definite symmetry with respect to the legs 
permutation (u-state and 5-state), and the triplet spectrum is 

E^ - 2A 2zJ PV' P^/4V'' 3Ji 3Ji K^ 

^„(,,„)-2A-2^„J(^-j ^-^j -_-±_^ + _. (25) 

Here M = 2m — 1/J is total mass of triplet excitation. 

The tunneling probability of the singlets from one leg to another is much smaller than that for the triplets. This is 
because the spinons with total spin zero cannot approach one each other: V^Cl) = (see discussion in the previous 

7 



section). Neglecting this tunneling we can say that the sing let g- and u-states are degenerate, and that the energies 
of the singlet states are given by ( [l2|) with substitution ( P4 ) . 

In conclusion to this section we would like to note that, for the ladder, the topological excitation with total spin 
or 1 is possible, as shown in Fig. 10. 



FIG. 10. Topological excitations with spins 1 or 0. 

However, the gap for this excitation is 2A and hence the threshold for creation of two topological excitations 
(minimum is two) is 4A. In the present paper we do not consider excitations with such high energies. 



IV. 2D ARRAY OF COUPLED MAJUMDAR-GHOSH CHAINS 



Another way of linking together spin chains is to form a two-dimensional array as shown in Fig. 11. 
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FIG. 11. The two-dimensional array. Each of the 'A' and 'B' chains represent spin chains similar to those previously 
considered. The short lines represent cross links with an exchange Jx. 

There are a total of N'^ cross links. The Hamiltonian is of the form 

X 

+ -J 2-^ \ x,y^x,y+l + ^■^^x,y^x,y+2) (26) 

y 

+ Jl- 2^ ^x,y^x,y 



The "A"-chains are aligned along the x-direction and the "B"-chains are aligned along the y-direction. First we want 
to demonstrate that the ground state of the system looks like that shown in Fig. 12. 
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FIG. 12. The two-dimensional array. The horizontal chains are the A-chains while the vertical chains are B-chains. The dots 
represent the cross- links, the same as those in Fig. 11. 



We stress that any dot in Fig. 12 denotes two spins 1/2 separated in the vertical direction (see Fig. 11) and coupled 
via J_L. The dimers aligned in a;-direction are built from A-chain spins and the diniers aligned in y-direction are built 
from B-chain spins. The question is: why is the ground state pattern like that in Fig. 13a (which is actually part of 
Fig. 12), but not like that in Fig. 13b? 
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(b) 
FIG. 13. Two arrangements of singlets in a two-dimensional array. 

At J_L = these configurations are degenerate. As was the case for the ladder (previous section), to calculate ground 
state energy correction due to J± it is convenient to use the localized triplet representation for virtual excitations. 
The energy of the localized triplet is J, and the amplitude of two triplets created at the two singlet bonds coupled via 
J± is J_l/4, see e.g. RefJlj. Therefore, the second order correction to the state Fig. 13a is described by the process 
|0) -^ |12) -^ |0) and is equal to 3(Jj_/4)^/(— 2J). Here |12) denotes the state with triplets excited at bonds 1 and 
2 (see Fig. 13a), the coefficient, 3, is number of triplet polarizations. However, a similar contribution exists for the 
configuration in Fig 13b and, therefore, the configurations Figl3a and Figl3b remain degenerate to second order in 
J± . To consider the next order we have to recall that the localized triplet can hop between two bonds coupled via J± . 
The hopping amplitude is also J_l/4, see e.g. Ref.ll3. Additional corrections for Fig. 13a arise due to the possibility of 
hopping around the square (closed loop trajectory), whereas, in Fig. 13b, there is no such possibility. The correction 
is due to the following sequences of virtual states (see Fig. 13a): 

(27) 
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The corresponding energy correction per one "dot" (one vertical J± cross-link) is then 

15(Jx/4)4 
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13a 



El3b 



J3 



(28) 



We have thus shown that the system prefers the formation of "squares" . However, even in this case there is a 
pattern shown in Fig. 14 which, in the considered order {Jj_), remains degenerate with that in Fig. 12. 
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FIG. 14. Another arrangement of singlets in a two-dimensional array. 



The degeneracy is lifted only in the J^ term. To see this we first must recall that the localized triplet can create an 
additional triplet at a nearby bond along the MG chain: 1'^^ -^ ^n<T*I+i,7 (*!« is the localized triplet creation operator 
at the bond n, (H is the triplet polarization). The amplitude of this process is jJiepaj, see Ref.S Keeping this in 
mind, one can see that the energy difference between Fig. 15a (building block of Fig. 12) and Fig. 15b (building 
block of Fig. 14) is due to closed loop trajectories of the type (for notations, see Fig. 15) 



|0) -^ |12) -^ |125> -^ |126) ^ |127) ^ |327) -> |427> -^ |42) -^ |0). 
These trajectories exist for the configuration in Fig. 15a but do not exist for Fig. 15b. 



(29) 





) 
2 6 

> 


1 f 








J I 


4 


7 
(a) 


■f. 








J e^ r 











r 




•1 









































U 



(b) 
FIG. 15. (a) Building block for Fig. 12a. (b) Building block for Fig 14. 

Calculation of the matrix element corresponding to the trajectory (g9|) and counting all the trajectories we find the 
energy difference. 
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(30) 



Thus, We have shown that the configuration presented in Fig. 12 is the true ground state of the Hamiltonian ([2(|). 
Note that the parameter of expansion in Eqs. (Pq ) and ( pOj) is J_l/4J. Note also that the coefficient in (|30|) is 
anomalously small meaning that the system is quite soft with respect to the dislocations shown in Fig. 14. 
The lowest energy excitation consists of two spinous, see Fig. 16. 
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FIG. 16. Two spinous along one of the A-chains in a two-dimensional array. 

Here the spinons are created on the upper A-chain and all other spins remains paired. According to (| 
tension of the confining string is 



T = 



15 (Ji/4)^ 



4 J3 

and hence the dynamics are described by Eqs. (|lO|), (^ and (|l^) with substitutions 

15 (Ji/4)4 
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We now intend to consider the spinons tunneling from one chain to another. This is simpler to do for triplet excitations 
and therefore we concentrate only on this case. Taking account of dl4) and (p2h, the wave functions of triplets built 
on the A- and B-chains are 
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where R = (ri + r2)/2, r = r2 — ri. Eqs. (33) do not make account of the tunnchng between the chains. The 
mechanism of the tunneUng is similar to that of the ladder (see Fig. 9) - th e sp inons approach one another to within 
one lattice space and then the perturbation, J^S^^Sb, in the Hamiltonian (Eq) swaps the triplet from the A-chain to 
the B-chain (and vice versa). Calculation of the tunneling amplitude gives 



Kx Ky A ^^^ Kx KyfJ^/A 

T = J± cos COS — 2. = 30 J COS cos — - 

2 2^^ 2 2 \ J 

The triplet dispersion is given by diagonalizing the 2x2 matrix 



where 



eA(K) = 2A - 2z„ J 



eB(K)=2A-2z„J 






3\'/VJ±\ 



'^ J3 + 2M 



8/3 



ijj 



J3 



e = t: ± 



£A - es 



2M 



Here M = 2m — 1/J is total mass of triplet excitation. Diagonalization of (Ba) gives 



(34) 

(35) 

(36) 
(37) 



(38) 



At very small momenta, K^/2Af -^ t oc J^, rotational invariance is restored and we have two distinct triplet 
excitations with a quadratic dispersion 



e « ±T — const 



2M 



±T. 



(39) 



V. CONCLUSIONS 

In the present work we have considered several models in the regime close to the deconfinement of spinons. These 
systems have complex spectra which consist of multiple singlet and triplet states. 

For the one dimensional Ji — J2 — (5 model with J2/J1 ~ 0.5 and (5 <C 1 we have calculated the splitting between the 
singlet and triplet states as well as the spectral weight of triplet excitations. The splitting is in a good agreement with 
available numerical data. On the other hand, the spectral weight is not consistent with results of series expansion. 
This may be explained by the very slow convergence of the series for the system which has a dense spectrum. It would 
be very interesting to study this problem using different numerical methods. 

For the Ji — J2 — J_l ladder (Jj_ —> 0) we found that the interval between levels is AE ex J^ . Each level has a 
fine structure: it consists of two triplet and two singlet states. The triplet-singlet splitting is oc J^. Splitting between 
the triplets is oc J^, which is due to tunneling between the legs. The typical size of the states under consideration is 

^ oc Jj_ . This is the radius of confinement. 

Similar results are obtained for the two dimensional array of coupled Ji — J2 chains (two dimensional Ji — J2 — J± 

model as J_i — > 0). The energy interval between the levels is AE oc J_l • Fine structure of the level consists of two 
singlets and two triplets. The triplet-singlet splitting is oc Jj_. Splitting between the triplets due to tunneling is oc J^. 

The confinement radius is ^ oc Jj^ .As does any lattice model, this model violates rotational invariance. However, 
in the low energy sector (K <C ( J_i/ J)^/^), rotational invariance is restored. 
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